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^ ■ Abstract 

o 

We consider quenches in non-conserved two-dimensional XY systems between 
0^ . any two temperatures below the Kosterlitz-Thouless transition. The evolving 



systems are defect free at coarse-grained scales, and can be exactly treated. 



i-rt ■ Correlations scale with a characteristic length L{t) oc t^" at late times. The 

O ■ autocorrelation decay exponent, A = (r/j + ?7/)/2, depends on both the initial 

and the final state of the quench through the respective decay exponents 

of equilibrium correlations, CEQ{r) ~ r~'^. We also discuss time-dependent 

quenches. 
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Phase-ordering kinetics is the study of the non-equihbrium process of equihbration after 
a rapid change of system parameters such as temperature or pressure ||1|. Typically, the 
system is quenched from a high-temperature disordered phase to a low-temperature ordered 
phase. The problem is challenging, because there are degenerate ground states competing 
to select the ordered phase, and because the evolution of the order-parameter is typically 
determined by a non-linear partial differential equation. Indeed, there have been only a 
handful of exact solutions in phase-ordering systems (see [||]): non-conserved scalar and XY 
[0(2)] systems in one dimension, and both non-conserved and conserved spherical models 
in general dimensions. Typically, phase-ordering systems have singular topological defects 
(such as domain walls, vortex lines, or point defects) seeded from the disordered initial 
conditions. While the structure of singular defects can be used to extract the growth laws 
of characteristic scales through Energy- Scaling arguments [§], their singular nature makes 
exact solutions of correlation functions difficult. In fact, the exact solutions mentioned 
above do not involve singular defects: one-dimensional XY systems have only non-singular 
topological textures |0], while spherical systems have no topological defects. It is natural, 
therefore, to consider other systems without topological defects. 

In this paper, we consider the coarse-grained two-dimensional (2D) XY model below 
the Kosterlitz-Thouless transition temperature Tkt, where there are no free vortices [Q. 
For non-conserved dynamics we exactly determine two-point correlations after a quench 
between any two temperatures at or below T^t- We find scaling solutions characterized by 
a single time- dependent length L(t) ~ t^^^, without the logarithmic factor [L ~ (t/lnt)^/^] 
associated with quenches to states with free vortices (i.e. from above T^t) 0)|5|. We also 
measure autocorrelations, for which the decay of the overlap with the initial conditions is 
characterized by an exponent A, through A{t) ~ L~^ 0^- For quenches to T = we 
find A = rii/2, where rji characterizes the initial asymptotic spatial correlations through 
C{r) ~ r~^\ This agrees with predictions of Bray et al. for quenches from a critical point 
to zero temperature [Q. For quenches between arbitrary temperatures below T^t "we find 
A = {rji +rif)/2, which depends on the initial anc? final states of the quench [rjf characterizes 



the asymptotic equilibrium correlations at the final state through CEQ{r) ~ r~^f]. We also 
discuss quenches with arbitrary temperature histories, and show that asymptotic correlations 
and autocorrelations are independent of the early temperature history. 

We consider over-damped, non-conserved, dissipative, "model A" dynamics with an equa- 
tion of motion 

dt$=-T5H/5$+^{^,t), (1) 

where 0(x) is a two-component order-parameter. The energy-functional, 

1 



H\d)]= / d^x 



^(V0)2 + V^(0) 



(2) 



has a potential with a symmetric global minimum at |0| = 1, such as V{(f)) = Vo(0^ ^ 1)^- 
The thermal noise, ^ = (^1,^2)5 is Gaussian distributed with zero mean, with correlations 
determined by the fiuctuation dissipation theorem: (^j(x, t)^j(x', t')) = 2TkBT6'^{'x—x')6(t — 
t')6ij. Below the Kosterlitz-Thouless transition for the 2D XY model 0] any vortices present 
will be bound in oppositely charged pairs. We only consider correlations at distances much 
larger than the characteristic pair size, so that we need only treat the renormalized spin- 
waves of the system. Effectively we work on the line of zero-fugacity renormalizat ion-group 
fixed points with < T < Tkt- In the limit Vq — ^ 00 the field, coarse-grained beyond 
the pair scale, has unit magnitude. We change to phase variables = e*^, and the energy- 
functional takes the form 

if = |/(V^)Vx, (3) 

where ps is the coarse-grained spin- wave stiffness of the system 0] . The equation of motion 
(0) simplifies to 

^k = -PsTk^e^ + a, (4) 

where we only keep the component of the thermal noise locally orthogonal to the order 
parameter 0, with (^k(^)^k'('^')) = '^^kBT6\t_-k_'6(t — t'). [We will adsorb V and ps into the 



time-scale, making [t] dimensionally equivalent to [l]"^ for the rest of the paper, except when 
we discuss quenches with general temperature history.] We have ignored the 27r periodicity 
of 6 in equation (^ because our system has no vortices after the coarse-graining, so that the 
phase can be taken to be everywhere continuous. 
The solution to the equation of motion (|p is 

^k(t) = ^k(0)e-''* + r rfte-^'(*-*)a(t), (5) 

^ u 

where the time t > is measured from the time of the quench. Since we start from an equili- 
brated state at or below Tkt, the initial phases are determined by the spin-wave Hamiltonian 
(0). The Fourier transformed phases are Gaussian distributed with a probability distribution 

P[R(0)}] (xexp|-5:^0k(O)0-k(O)|. (6) 

We use rji and rjj to describe the initial and final quench states, respectively, of our system, 
with 1]{T) = ksT /2Tips{T). They describe the decay of equilibrium correlations and so are 
directly measurable in experiments. The "temperature" T always indicates the combination 
of system parameters (temperature, pressure, composition, . . . ) that determines ps and rj. 

The phase-phase correlations at general times after a quench at t = from a temperature 
Tj to a temperature Tj, both at or below Tkt 0]' ^^^ then straight-forward to calculate 
from d^): 

In particular, the initial correlations are given by (^k(O)^-k(O)) = 2'Krii/k'^. The phase- 
difference correlations then follow, 

i?(r,t,t') = ([^(x,t)-^(x + r,t')F 



. . 2 [(ek(t)^-k(t)) + (^k(t')^-k(t')) - 2 cos (k • r) (^k(t)^-k(t'))] 
= BEQir,t,t') + BNEQir,t,t'). (8) 

Beq and Bneq are the equilibrium and non-equilibrium correlations at the final tempera- 
ture, given by the terms proportional to rjf and rji — rjf, respectively, in (j^): 



BEQ{r,t,t') = Vf{l + \n{ryAal) + E,[ry4{al+\t-t'\)]}, (9) 

and 



^2 



\4^{al + 2t)ial + 2t')) J 

where Ei(x) = J^dye~'^/y for x > and 7 ~ 0.577 is Euler's constant. [We have used 
a soft ultra-violet cutoff, through a factor exp(— agA;^) in the integrand of each k-integral, 
where oq is of order the lattice spacing.] Both Beq and Bneq are manifestly symmetric 
under interchange of t and t'. Beq only depends on the magnitude of the time-difference 
\t — t'\, as expected for an equilibrium correlation. Conversely, B^eq has a scaling form at 
late times, 

BNEQir, t, t') = F{r/Vi, r/Vf), t, t' > a^ (11) 

where F{x) is a time- independent scaling function. 

Order-parameter correlations follow directly from the phase correlations since the phase 
variables are Gaussian distributed at all times — due to the linear evolution equation (H), 
and the Gaussian nature of the noise and the initial conditions. The general two-point 
two-time correlation function is 



C(r,t,t') = (</>(x,t)-</)(x + r,t') 

= (cos([e(x,t)-0(x + r,t')])), 

= exp{-E(r,t,t')/2} 

= CEQ{r,\t-t'\)C^EQ{r,t,t'), (12) 

where Ceq^t, \t-t'\) = exp{-BEQ{r,t,t')/2} and CNEQ{r,t,t') = exp {-B]yEQ{r,t,t')/2}. 
We see that C{r, t, t') has a product scaling form, with a growing length-scale L oc t^'"^ 
characterizing the non-equilibrium factor Cneq through Eq. (pAl). For quenches to Tf > 
0, the equilibrium factor Ceq has a non-trivial form equal to the equilibrium correlation 
function of the critical point at the final temperature. This is a generalization of the scaling 



expected for a quench to within an ordered phase, dominated by a T = fixed point, where 
Ceq{i^) = (101)^- We do not expect the product form seen in Eq. (|T2|) to hold in general 
0{n) systems for quenches to critical points; the product form holds for this XY system due 
to the Gaussian distribution of the phase variables. 

It is interesting to contrast these scaling results to the scaling prediction L ~ (t/lnt)^/^ 
0, for quenches with free vortices, observed in simulations ^ (note also [|l^). Quenches 
from below Tkt o^^e consistent with previous Energy-Scaling predictions , which are based 
on the observed late-time defect structure for quenches to T = 0. For non-conserved 
quenches without topological defects, a growth law L(t) ~ t^/^ is predicted for systems 
which scale - in agreement with the results of this paper. For quenches to T > 0, the 
Energy- Scaling approach does not directly apply. However, it is reassuring that scaling and 
the same growth law is observed in the non-equilibrium correlations Cneq- 

We determine the asymptotic correlations by using the asymptotics of Ei(x): 

{—7 — Inx, X ^ 1, 
(13) 
e~^ /x, X ^ 1. 

These determine the asymptotic equilibrium correlations, 

CsQ(r,0)~(r/ao)-''^ r » ao, (14) 

which reproduces the standard result 0]. The full equal-time correlations after the quench 
have the asymptotic behavior 

C{r,t,t)c^l (15) 

[(r/ao)-"/, rVt<l, 

where r, y/t ^ oq. They have the same spatial dependence as the equilibrium correlations 
(with the final temperature determining the correlations at short distances and the initial 
conditions determining them at long distances), but have an additional amplification factor 
at long scales. For rji = rjj we recover the equilibrium correlations (p^, as expected. 

The autocorrelation function is given by the r = correlations with the initial conditions. 
From equations (|]) to ([T3|), 



/ - \ -(Vi+Vf)/'i: 

A(t) = C(0,t,0)~(^j , (16) 

where again we take t ^ a^. Using L ~ t^/^, and A{t) ~ L^'^ |^, we determine the exponent 
^ = iVi ~^ ^/)/2 describing the decay of autocorrelations after a quench. Interestingly, A 
does not change if the quench is reversed: from rjf to T]i. For the decay of equilibrium 
autocorrelations, we set rjf = rjf = r] and find \eq = f]. 

Bray et al. have considered quenches from long-range correlated initial conditions 
to Tj = {rjf = 0). They predict that if initial long-range correlations are described by 
C(r, 0) ~ r"-'^, and if a < a^, then C(r, t) ~ {L/ry-'' for r > L and A{t) ~ t('^-'^)/4 for large 
t. For quenches from a critical point, cr = 2 — 77. Our results for quenches to T = agree 
(with d = 2): from Eq. (pIS]) we have C ~ (L/r)'^^ with L ~ t^/^, while from Eq. (^) we have 
A(t) ~ t^''»/^. The initial correlations are relevant for quenches from all temperatures at or 
below Tkt, so that Cc < 2 — rj{TKT) = 7/4. Because Xsr, the value for short-range correlated 
initial conditions, is never less than A for long-range correlated initial conditions |0 , we have 
^SR > ^(Tkt) = 1/8 for a quench to Tj = 0. This lower bound is much smaller than the 
values measured by Pargellis et al. f^ . Theoretical treatments of quenches to critical points 



14 have previously only treated the case of a quench with uncorrelated initial conditions, 
where A is independent of the details of the initial conditions. The methods of reference 
0] can be extended to general systems, in d-dimensions, with a quench to a general critical 
point Tf (characterized by an equilibrium correlation function C{r) ~ r^^'^^'^'^'^f^), from a 
state with power-law spatial correlations (in the same order-parameter field) decaying as 
^-(d-2+r},)^ We find A = rf — 2 + (r^j + ?7/)/2 |13[, provided that the long-range correlations in 
the initial state are relevant, which requires that the exponent Xsr for a quench to Tf from 
initial conditions with short-range correlations satisfy Xsr < d+ {rjf — rii)/2. For quenches 
between two critical points in the 2D XY model, our result A = (77^ + ^/)/2 agrees with the 
general result for d = 2, and depends on both the intitial and the final temperatures of the 
quench. 

Arbitrary temperature histories of the quench can be treated in a similar manner to 



our discussion so far. We assume that the thermal bath has a well-defined time-dependent 
temperature T = T(t), i.e. that microscopic time-scales are much faster than the phase- 
ordering time-scales. Then the renormalized spin-wave stiffness, Ps{T), and the equilibrium 
decay exponent, ri(T), will both be time- dependent through their temperature dependence. 
Eq. (|I|) to (^) will be unchanged, while Eq. (^ will change to 

^k(t) = 0k(O)e-'='p(*) + /* rfte-'='[p(*)-P(*)]a(tl, (17) 

Jo 

where p(t) = /q ps(t)dt. This leads to phase correlations 

(^k(t)^-k(t')) = ^e-'^'lPW+P^*')] + 2kB Z""^*'*'^ rfte-'='[PW+P(*')-2p®lT(t). (18) 

k'^ Jo 

This can be used in Eqs. (|]) and (|12D to determine correlations under an arbitrary tem- 
perature history. It is straight-forward to show that temperature changes before a time tM 
do not affect correlations with t,t' ^ tM or autocorrelations A{t) for t ^ Im- This is best 
illustrated by considering autocorrelations in the low-fugacity limit, where ps is temperature 
independent. With some work, we have 5(0, 0, t) = (77^/2) In (t/al) + Jq dir]{i)/{al + 2t- 2i). 
[We have adsorbed the constant ps into ri{t).] If we restrict the time-dependence of T to 
times before tju -C t, then _B(0,0,t) = In (t/al){rii + ?7/)/2 + 0{tM/t) and we recover our 
previous result. The same approach demonstrates that late time correlations and autocorre- 
lations are insensitive to the temperature history before a time Im <^ t- Exact correlations 
for simple quench histories can easily be calculated, and may be useful for comparing to 
experiments, but do not appear qualitatively different than instantaneous quenches. 

The ease of calculation of this model is fortuitous. In part this is because, for vector 
systems, the gradient term in equation (0) dominates over the potential term p[. Since the 
potential term is not needed to set a core-scale for singular defects (such as free vortices), 
the hard-spin limit may be taken without complications. The 2D XY model is special below 
Tkt because the vortices are tightly bound and only serve to renormalize the spin-wave 
stiffness ps of the effective hard-spin Hamiltonian. The resulting Gaussian nature of the 
phase variables greatly simplifies the analysis. [Note that the hard-spin limit for conserved 



"model B" dynamics leads to a much more complicated evolution equation than (^), and 
the phases will not be Gaussian distributed.] Of course our results only apply for distances 



large with respect to the vortex pair size in the two equilibrium phases at Tj and Tf p^|. 

The specially prepared nematic system of Pargellis et al. 0, developed to exhibit 2D XY 
behavior, should exhibit the behavior described by this paper. The experimental procedure 
will be easier since late-time correlations and autocorrelations do not depend on the early 
stages of the quench. [Experimental analysis will not be able to rely, however, on the 
characteristic schlieren patterns of free vortices, which will be absent.] Most other planar 
nematic liquid crystal systems, while they exhibit an XY-like growth law when free vortices 
are present 0, will not have their vortex- free quenches described by this paper. This is 
because the spin- waves of an unconstrained nematic order parameter are not of the 0(2) 
variety - the spins are not characterized by a single angle. 

We have calculated correlations for quenches in the 2D XY model between any two 
phases at or below the Kosterlitz-Thouless temperature. The non-equilibrium part of the 
correlations scales with a characteristic length-scale L{t) ~ t^^^. This growth law differs by 
a logarithmic factor from the growth laws expected for quenches involving free vortices. The 
autocorrelation decay exponent depends on both the initial and the final state of the quench, 
A = {rii + T]f)/2. The asymptotic autocorrelations and equal-time correlations do not depend 
on the early temperature history of the quench, but only on the initial conditions (through 
rji) and on the temperature at late times (through rjf). 

We thank T. Blum and D. A. Huse for discussions, and the Isaac Newton Institute for 
hospitality. 
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